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ABSTRACT 


The  classical  Smith  Predictor  method  for  single  variable 
systems  with  time  delays  is  extended  to  a  class  of  linear  multi  - 
variable  systems.  Derivations  of  the  multivariable  Smith  Predictor 
are  presented  for  both  continuous  time  and  discrete  time  systems 
which  contain  time  delays  in  the  control  variables  and/or  output 
variables.  As  in  the  classical  method,  use  of  the  multivariable 
Smith  Predictor  eliminates  the  time  delays  from  the  characteristic 
equation  of  the  closed-loop  system.  The  multivariable  Smith 
Predictor  is  applied  to  a  double  effect  evaporator  pilot  plant  in 
both  simulation  and  experimental  runs  and  the  effect  of  modelling 
errors  on  the  predictor  response  is  also  examined  in  the  simulation 
runs . 

Finally,  algorithms  are  derived  for  more  complex  cases 
such  as  systems  with  only  some  delayed  control  variables  or  with  time 
varying  or  inaccurately  determined  delays. 

The  simulated  and  experimental  results  demonstrate  the 
ability  of  the  multivariable  Smith  Predictor  method  to  efficiently 
handle  delays  in  the  control  and/or  output  variables  while  main¬ 
taining  freedom  in  the  design  of  the  feedback  control  matrix. 
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CHAPTER  ONE 
INTRODUCTION 

The  existence  of  time  delays  is  of  common  occurrence  in  many 
process  control  systems.  In  fact,  classical  methods  for  process 
modelling  often  include  a  time  delay  in  the  system  transfer  function. 

Time  delays  often  occur  due  to  transportation  lines  for 
fluids,  pneumatic  instrument  signals,  etc.  or  due  to  the  time  delays 
associated  with  measurements  (e.g.  time  required  to  analyze  chemical 
compos i ti on) . 

The  detrimental  effects  of  time  delays  on  system  stability 
and  control  are  well-known  to  both  control  system  designers  and  per¬ 
sonnel  responsible  for  plant  operation.  From  the  classical  viewpoint, 
the  phase  lag  introduced  by  a  time  delay  tends  to  reduce  system 
stability  and  make  satisfactory  control  more  difficult  to  achieve. 
Furthermore,  time  delays  greatly  complicate  control  system  design  for 
multivariable  systems  since  many  design  approaches  are  either  not 
applicable  to  time  delay  systems  or  the  resulting  control  system  may 
be  unduly  complicated  and  difficult  to  apply. 

1 . 1  Literature  Survey 

Several  methods  for  the  control  of  multivariable  systems  with 
time  delays  have  been  proposed  in  the  literature.  For  example,  opti¬ 
mal  control  design  methods  have  been  extended  to  systems  with  time 
delays  [1-4],  but  except  for  special  cases,  the  derivation  of  the 
control  law  is  often  difficult,  relies  heavily  on  complicated  control 
algorithms  and  consequently,  actual  implementation  of  these  optimal 
control  algorithms  in  process  control  systems  is  seldom  attractive. 
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On  the  other  hand,  various  methods  for  single  variable  control 
systems  have  been  reported  in  the  literature  that  utilize  a  predictive 
approach  for  compensation  of  the  time  delays,  but  it  appears  that  these 
methods  have  been  applied  to  actual  processes  in  only  very  few  cases 
[5-12].  These  predictive  methods  involve  the  use  of  a  time  delay  in  a 
feedback  loop  around  the  controller  since  a  controller  which  uses  the 
information  that  the  process  has  a  delay  should  be  able  to  make  more 
intelligent  corrections  than  a  controller  that  receives  only  the  delayed 
error  signal . 

An  efficient  continuous  time  predictor,  (the  Smith  Linear 
Predictor),  was  derived  by  Smith  [5,6],  for  single  variable  systems. 
Applications  of  this  predictor  scheme  and  comparative  studies  with 
classical  I,  PI  and  PID  controllers  were  reported  in  the  literature 
[7-10],  for  continuous  and  sampled  data  systems.  The  conclusions 
resulting  from  these  studies  do  not  agree  in  all  the  cases  especially 
as  far  as  the  effectiveness  of  the  Smith  Linear  Predictor  for  regulatory 
control  is  concerned.  Shinskey  [7],  and  Nielsen  [8]  reported  that  the 
Linear  Predictor  and  the  similar  "complementary  feedback  technique"  are 
of  questionable  value  for  regulatory  control  of  time  delay  processes 
while  Buckley  [9,10]  and  Lupfer  and  Oglesby  [11,12]  applied  the  Smith 
Linear  Predictor  on  actual  processes  with  excellent  results  for  regula¬ 
tory  control.  Single  variable,  sampled  data  controllers  utilizing  the 
Smith  Predictor  approach  have  been  designed  and  simulated  results  have 
been  reported  in  the  literature  [13,14].  Moore  [15]  derived  a  digital 
algorithm  for  single  variable  discrete-time  systems  with  time  delays 
that  is  conceptually  similar  to  the  Smith  Linear  Predictor,  the  differ¬ 
ence  being  mainly  a  more  extensive  use  of  the  mathematical  model  of  the 
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process  than  in  the  Smith  Predictor  scheme.  Results  for  actual  pro¬ 
cesses  with  this  predictor  indicate  a  satisfactory  regulatory  control 
and  better  performance  than  the  classical  methods  of  I,  PI  and  PID 
control  [15]. 

Although  many  studies  have  considered  predictors  for  single 
variable  systems,  no  suitable  predictor  has  been  reported  for  multi - 
variable  systems  with  time  delays.  The  analysis  required  presents 
more  difficulties  and  complexities  than  in  the  single  variable  case 
since,  typically,  combinations  of  delayed  and  undelayed  process 
variables  and/or  measurements  occur  in  multivariable  control  systems. 

Jacobson  [16]  proposed  an  extension  of  the  digital  predictor 
algorithm  derived  by  Moore  [15]  to  compensate  for  the  measurement 
delay  introduced  by  the  digital  interface  in  computer  control  systems. 

He  applied  this  scheme  to  a  double  effect  evaporator  pilot  plant  at 
the  University  of  Alberta  with  satisfactory  results;  however,  no  general 
method  has  been  reported  for  the  more  general  case  of  time  delays  in 
multivariable  control  systems. 

1 .2  Objectives  of  the  Study 

This  study  is  concerned  with  the  development  of  a  suitable 
predictor  scheme  for  the  control  of  multivariable  systems  that  contain 
time  delays.  The  Smith  Linear  Predictor  was  extended  to  the  case  of 
multivariable  systems.  Since  the  Smith  Linear  Predictor  permits  design 
of  the  control  by  using  methods  suitable  for  undelayed  systems,  this 
method  was  considered  particularly  suitable  for  multivariable  systems. 

The  derivation  of  the  control  algorithm  had  to  conform  to  the 
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(i )  Satisfactory  control  for  linear  multivariable  systems  with  time 
delays  in  the  manipulated  and/or  the  measurement  variables. 

(ii)  Ability  to  handle  the  cases  when  only  some  of  the  control  vari¬ 
ables  or  some  of  the  measurement  variables  are  delayed. 

(iii)  Easy  implementation  of  the  control  algorithm  on  a  digital  control 
computer. 

1.3  Structure  of  the  Thesis 

The  thesis  is  concerned  with  the  derivation  and  evaluation  of 
a  predictor  algorithm  for  several  classes  of  multivariable  systems  with 
time  delays. 

In  Chapter  Two,  the  theoretical  derivation  is  presented  for 
certain  continuous  and  discrete  time  systems  together  with  a  theoretical 
investigation  of  more  complex  cases  and  techniques  for  situations 
encountered  often  in  practice.  In  Chapter  Three,  simulated  results  of 
the  multivariable  predictor  are  given  for  the  double  effect  evaporator 
model  and  several  combinations  of  measurement  and  process  time  delays. 

In  Chapter  Four,  the  results  of  an  experimental  application  of  the 
multivariable  predictor  on  the  evaporator  pilot  plant,  using  the  IBM  1800 
process  control  computer  are  given  and  in  Chapter  Five,  the  overall 
conclusions  are  given  for  the  multivariable  predictor  control  scheme 


developed  in  this  thesis. 
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CHAPTER  TWO 

THEORETICAL  DEVELOPMENT 

2.1  Introduction 

In  1957,  O.J.M.  Smith  (Smith  1957,  1959)  proposed  a  control 
technique  for  single  variable  control  systems  which  contain  time  delays. 
This  technique,  which  became  known  as  the  Smith  Predictor  (or  Smith 
Linear  Predictor),  is  illustrated  in  Figures  2.1  and  2.2.  The  chief 
advantage  of  the  Smith  Predictor  method  is  that  time  delays  are 
eliminated  from  the  characteristic  equation  of  the  closed-loop  system. 
This  is  achieved  by  including  a  mathematical  model  of  the  process  in 
the  feedback  loop  around  the  controller. 

The  output  of  the  predictor  block  in  Figure  2.2  is  the  dif¬ 
ference  between  two  model  responses:  the  response  of  the  system  with¬ 
out  the  time  delay  minus  the  response  of  the  system  with  the  time 
delay.  If  the  process  models  were  perfect,  then  the  actual  process 
response,  y(s),  would  be  cancelled  by  the  model  response  and  the 
control  action  would  be  based  on  the  response  of  the  model  without  time 
delay.  For  the  control  system  in  Figure  2.2,  the  closed-loop  transfer 
function  for  load  changes  is 


(s)  r  ,  ,  Gc(s)Gp(s)GL(s)e 
(it  =  VS)  ‘  1  +  GWs)G  (s)H 


x(s 

d 


-as 

JH Is) 


(2.1) 


c'  '  p 

The  characteristic  equation  for  the  closed-loop  system  in  Figure  2.2  is 
gi ven  by 


1  +  Gc(s)Gp(s)H(s)  =  0 


(2.2) 
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Figure  2.1  Single  Variable  Control  System  with  a  Time  Delay  in  the 
Process  Transfer  Function 


Predictor 


Figure  2.2  Smith  Predictor  for  the  Feedback  Control  of  Figure  2.1 
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which  is  also  the  characteristic  equation  for  the  system  in  Figure 

2.1  when  the  time  delay  is  zero.  Thus  the  Smith  Predictor  has  success¬ 
fully  eliminated  the  time  delay  from  the  characteristic  equation  and 
the  design  of  a  suitable  control  matrix  can  be  done  with  any  method 
appropriate  for  undelayed  systems. 

2.2  Multivariable  Continuous  Time  Systems 

Consider  the  following  linear,  stationary,  state-space 
model  with  time  delays  in  the  output  variables  and  the  control  vari¬ 
ables: 

x ( t )  =  A  x (t )  +  B  u(t  -  a)  +  D  d(t)  (2.3) 

y(t)  =  C,  x(t)  +  C2  x(t  -  b)  (2.4) 

where:  x(t)  =  state  vector  of  dimension  n, 

u(t)  =  control  vector  of  dimension  m, 
d(t)  =  disturbance  vector  of  dimension  p, 
y (t)  =  output  vector  of  dimension  r, 
a,b  =  constant  time  delays 

A>B,C-|,C2  and  D  are  constant,  real  matrices  of  appropriate 
dimensi ons . 

The  time  delays  in  this  state-space  model  can  be  given  the 
following  physical  interpretation.  Time  delay  a,  is  associated  with 
the  calculation  and  implementation  of  control,  that  is,  a  time  delay 
in  all  control  variables.  The  more  general  case  of  some  control  vari¬ 
ables  delayed  and  some  not  is  theoretically  examined  in  §  2.4.  In 
many  physical  systems,  time  delays  are  also  associated  with  the  measure¬ 
ment  of  certain  state  variables.  A  notable  example  in  process  control 


I 

z  •  ■  7  ?ri:  ,  v:.“ 
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systems  is  chemical  composition.  This  variable  is,  in  general,  diffi¬ 
cult  to  measure  and  often  requires  a  period  of  time  to  carry  out  the 
analysis,  i.e.  a  time  delay.  The  inclusion  of  time  delay,  b,  in 
Equation  (2.4)  provides  a  general  model  for  systems  in  which  some 
state  variables  can  be  measured  instantaneously  but  time  delays  are 
involved  in  the  measurement  of  other  state  variables.  Apparently,  the 
output  equation  in  Equation  (2.4)  has  not  been  widely  used  in  previous 
investigations  of  time  delay  systems  despite  its  obvious  practical 
importance. 

2.2.1  Time  Delay  in  Output  Variables 

If  no  time  delay  is  present  in  the  control  variables  (i.e. 
a  =  0  in  Equation  (2.3)),  the  state-space  model  in  Equations  (2.3)  and 
(2.4)  reduces  to 


x(t)  - 

A  x(t)  +  B  u(t)  +  D  d(t) 

(2.5) 

J£(t)  - 

C1  x(t)  +  C2  x(t  -  b) 

(2.6) 

Assuming  zero  initial  conditions  and  taking  the  Laplace  transform  of 
Equations  (2.5)  and  (2.6)  gives,  after  rearrangement, 


x(s)  -  G  (s)  u(s)  +  Gj^(s)  d(s) 

(2.7) 

y(s)  =  C,  x(s)  +  C2  e"bs  x(s) 

(2.8) 

where  j  =  nxn  identity  matrix 

Gp(s)  =  (si  -  A)~^B  =  process  transfer  function  matrix 

Gl(s)  h  (si  -  A)_1D  =  load  transfer  function  matrix 
and  it  is  implicitly  assumed  throughout  this  thesis  that  matrix  inverses 


exist. 
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Suppose  an  output  feedback  control  law  of  the  form 

u(s)  =  -  Gc(s)  y(s)  (2.9) 

is  assumed  where  G  (s)  is  the  matrix  of  feedback  controller  transfer 

=c 

functions.  Then  combining  Equations  (2.7)  -  (2.9)  gives  the  following 
characteristic  equation: 

I  In  +  Gp(s)Gc(s) {£i  +  C2e_bs)  I  =  0  (2.10) 

where  I  is  an  n  x  n  identity  matrix  and  the  symbol  |  |  denotes  the 
determinant. 

Next,  it  will  be  demonstrated  that  the  Smith  Predictor  method 
can  be  used  to  eliminate  the  time  delay  from  the  characteristic  equation. 
Consider  the  block  diagram  shown  in  Figure  2.3.  As  for  the  single 
variable  system  in  §  2.1,  the  feedback  loop  around  the  controller  con¬ 
tains  a  mathematical  model  of  the  process,  both  with  and  without  time 
delays.  From  the  block  diagram  in  Figure  2.3,  it  follows  that  for 
r(s)  =  0, 

u(s)  =  -  Gc  *(s)  -  Gc  C2  Gp(l  -  e'bs)  u(s)  (2.11) 

where  Gc  and  Gp  denote  Gc(s)  and  Gp(s),  respectively.  Rearranging 
Equation  (2.11)  gives  the  following  control  algorithm 

u(s)  =  §c[^(s)'?25pe"bS  y(s)]  (2J2) 

Combining  Equations  (2.11),  (2.7)  and  (2.8)  gives 

x(s)  =  GLd(s)  -  Gp^C^ts)  -  GpGcC2(x(s)-GLd(s)) 

-  =p=c=2e"bS=L-(s) 


(2.13) 


j  iFnr:-H3j so 

(f-r.s)  (t)y  (ed*s  -  r)  §  so  oa  -  Uh  ,3  -  *  («)«  '  9 


j  •  b  (v  s)  «  rr.- )  e  f j. i 

| 

((*:  *  <•>**  ‘  <*>« 
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Figure  2.3  Smith  Predictor  for  a  Multivariable  System  with  a  Time  Delay 
in  the  Output  Variables 
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or  rearranging  gives 

x(s)  =  Cln+GpGc (C-|+C2 ) ]  1  [In+GpGcC2(l-e  bs)]GLd(s)  (2.14) 

From  Equation  (2.14)  it  is  apparent  that  the  characteristic  equation 
for  the  closed  loop  system  is 

I  In  +  Mc(=Cl  +  =C2)!  =  0  (2J5) 

which  is  also  the  characteristic  equation  for  the  system  without  time 

delays  (i.e.  b  =  0  in  Equation  (2.6)).  Thus  the  multivariable  Smith 

Predictor  has  eliminated  the  time  delay  from  the  characteristic 

equation.  The  design  of  the  controller  transfer  function  matrix 

G  (s)  can  then  proceed  using  design  techniques  developed  for  systems 
“  c 

without  time  delays. 

2.2.2  Time  Delays  in  Both  Control  Variables  and  Output  Variables^ 

The  case  is  considered  where  constant  time  delays  occur  in 
both  the  control  variables  and  the  measured  outputs.  The  system  of 
interest  is  given  in  Equations  (2.3)  and  (2.4).  Assuming  zero  initial 
conditions  and  taking  the  Laplace  transform  of  Equations  (2.3)  and 
(2.4)  gives 


x(s)  =  Gp  e  as  u(s)  +  Gl  d(s) 

(2.16) 

y(s)  =  G-j  x(s)  +  C2  e  x(s) 

(2.17) 

The  multivariable  Smith  Predictor  for  this  system  consists  of 
two  feedback  loops  around  the  controller  as  shown  in  Figure  2.4.  From 
Figure  2.4  it  follows  that  for  r(s)  =  0, 


u(s)  =  -^(s)  -  §C%+S2^=p-^ 


G  C-i  G  e 
=c=l=p 


-as 


u(s)  +  GcC2Gpe'(a+b)su(s)  (2.18) 
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Figure  2.4  Smith  Predictor  for  a  Multivariable  System  with  Time  Delays  in  Both 
Control  and  Output  Variables 
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or  rearranging  gives  the  following  control  algorithm 


u(s)  =  -[In+Gct^+C^Gp]'1  Gc[y(s)-(C1-C2  e 


(2.19) 


Combining  Equations  (2,16),  (2.17)  and  (2.18)  gives  the  result 


From  Equation  (2.20)  it  follows  that  the  characteristic  equation  for  the 
closed  loop  system  is 


(2.21) 


In  +  SpSc(£l^z)l  =  0 


=p=cv=l  =2 


Thus,  the  Smith  Predictor  is  again  successful  in  eliminating  time 
delays  from  the  characteristic  equation.  Furthermore,  the  resulting 
characteristic  equation  in  Equation  (2.21)  can  easily  be  shown  to  be 
the  characteristic  equation  for  the  system  without  time  delays  (i.e. 
the  system  in  Equations  (2.3)  and  (2.4)  with  a  =  b  =  0). 

2.3  Multivariable  Discrete-Time  Systems 


A  multivariable  Smith  Predictor  can  be  developed  in  an 


analogous  fashion  for  discrete-time  systems  which  contain  time  delays. 
Consider  the  following  stationary  state-space  model 


x(n+l)  =  $  x(n)  +  e  u(n-a)  +  a  d(n) 


(2.22) 


y(n )  =  C-|  x(n)  +  C2  x(n-b) 


(2.23) 


where  x(n)  =  state  vector  of  dimension  n. 


u(n)  =  manipulated  vector  of  dimension  m, 
d(n)  =  disturbance  vector  of  dimension  p, 


( b  !  •  3  '  (f.  >  ’  f  +  . . 
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y(n)  =  output  vector  of  dimension  r, 

a,b  =  constant  time  delays  which  are  integers  (i.e.  positive 

integer  multiples  of  the  sampling  period,  T) 

4,e,A,C.j  and  C2  are  constant  matrices  of  appropriate  dimen¬ 
sions. 

As  in  the  continuous  time  case,  a  multivariable  Smith 
Predictor  algorithm  can  be  developed  by  considering  feedback  loops 
around  the  controller.  Suppose  a  proportional  feedback  control  law 
is  assumed  of  the  form 

u (n )  =  -  Kc  y(n)  (2.24) 

where  K  is  a  constant  m  x  r  matrix.  Then  a  suitable  algorithm  for  the 
=c 

Smith  Predictor  is 


u(n)  =  -  Kc  y(n)  -  Kcp(n)  (2.25) 

p(n)  =  C-jP-j  (n)  +  C2P2(n)  (2.26) 

where 

Ej  (n)  =  <£  gj(n-l)  +  §  u(n-l)  -  9  u(n-a-l)for  n  1  1  (2.27) 

g2(n)  =  4>02(n“l)  +  6  u (n- 1 )  -  9  u(n-a-b-l)  for  n  ^  1  (2.28) 


and 

p-j  (0)  =  p2(0)  =  0 


This  formulation  of  the  predictor  algorithm  eliminates  the  time  delays 
from  the  characteristic  equation  of  the  closed-loop  system  as  it  is 
shown  below.  Taking  the  z-transform  of  Equations  (2.22),  (2.23)  and 


(2.25)  -  (2.28)  gives 


*  o  it  "  -  J  f 
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z  X(z) 

=  Mez"a  U(z)  +  MaD(z) 

(2.29) 

Y(z) 

=  C^tz)  +  C2z"°X(z) 

(2.30) 

y(z) 

=  KcY(z)  -  KcP(z) 

(2.31) 

P(z) 

=  c1P1(i)  +  £2P2(z) 

(2.32) 

P-|(z) 

=  z'J^fz)  +  eU(z)  -  ez"aU(z)] 

(2.33) 

P2(z) 

=  z'[|P2(z)  +  eu(z)  -  z‘(a+b)8U(z)] 

(2.34) 

where,  the  notation,  X(z)  =  z- transform  of  x(n),  is  used.  Equations 
(2.29)  -  (2.34)  can  be  combined  and  rearranged  to  give  the  following 

expression : 

x(z)  ^zl^Me^C^+CgirHl^MeK^O-z"8)  z'1 

+  Me KcC2 (1 “ Z_ ^a+b ^ ) z- ' ]  MfiD(z)  (2.35) 

where  M(z)  is  defined  as 

M(z)  =  (  In  -  z-1!)’1  (2-36) 

Equation  (2.35)  implies  that  the  characteristic  equation  for  the  system 
with  the  Smith  Predictor  is 

|  z  !n  +  M9Kc(C)  +  C2)  I  =  o  (2.37) 

It  can  easily  be  shown  that  Equation  (2.37)  is  also  the  characteristic 
equation  for  the  original  system  in  Equations  (2.22)  and  (2.23)  when 
the  time  delays  are  zero  (i.e.  a  =  b  =  0).  Thus  the  Smith  Predictor 
algorithm  in  Equations  (2.25)  -  (2.28)  can  be  used  to  eliminate  time 
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delays  from  the  characteristic  equation.  The  design  of  the  feedback 
controller  matrix  K  can  then  proceed  using  the  wide  variety  of  design 
techniques  developed  for  systems  without  time  delays.  The  above 
analysis  could  also  be  extended  to  other  types  of  multivariable  control 
techniques  such  as  proportional  plus  integral  control. 

2.4  Application  of  Predictor  Algorithm  to  Systems  With  Delay  Only  in 

Some  Manipulated  Variables 
2.4.1  Introduction 

The  case  of  systems  with  some  control  variables  delayed  and 
some  not  can  be  considered  a  generalization  of  the  previously  considered 
case  of  a  delay  in  all  of  the  control  variables.  From  a  practical  point 
of  view  this  situation  is  more  commonly  found  in  practice  where  usually 
it  is  expected  that  the  delays  in  the  various  control  variables  will 
be  different  or  zero  (different  length  of  pneumatic  transmission  lines, 

etc. ) . 

/ 

Assume  that  a  system  can  be  described  by  the  following  state- 
space  model : 


x(t)  =  gx(t)  +  B-j  u^t)  +  B2  u2(t-a) 


(2.38) 


y(t)  =  x(t) 


(for  simplicity) 


(2.39) 


Here,  the  control  vector  has  been  partitioned  as 


u-j  ( t) 


u(t) 


(2.40) 


u2(t) 


where  u-j  denotes  the  i  undelayed  control  variables  and  u2  the  (m-£) 


. 

. 


- 
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control  variables  which  are  delayed.  An  equivalent  representation  of 
the  system  in  Equations  (2.38)  and  (2.39)  is  given  by 


x  ( t )  =  Ax(t)  +  B-jU(t)  +  §2u(t-a) 
y(t)  =  x(t) 


where  B-j  and  are  defined  by 


=B1 


(2.41) 

(2.42) 


can 

1 _ 

r~ 

O  II 

*  _ 

(2.43) 


=2 


=  2 


0 


(2.44) 


where  B^  and  §2  are  n  x  m  matrices  and  the  notation  is  the  same  as  in 
§  2.2  for  the  other  matrices  and  vectors. 

Taking  the  Laplace  transform  of  Equations  (2.41)  and  (2.42) 


gi  ves 


where 


x($)  =  G^sjuts)  +  G2(s)e'asu(s)  +  G2(s)d(s) 
y (s )  =  x(s) 


(2.45) 

(2.46) 


-1 


G-,  (s)  5  (si  -  A)  B-j 


-1 


G2(s )  -  (si  -  A)  B2 


(2.47) 

(2.48) 


In  applying  the  predictor  scheme  on  a  system  described  by 
Equations  (2.45)  and  (2.46)  a  difficulty  is  encountered  due  to  the 
interaction  between  the  delayed  and  undelayed  parts  of  the  system. 
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The  search  for  a  suitable  predictor  scheme  for  this  system 
was  directed  towards  finding  an  algorithm  that  will  in  effect  decouple 
the  system  in  the  sense  that  any  predictive  action  on  the  delayed 
manipulated  variables  will  not  affect  the  undelayed  portion  of  the 
system. 

2.4.2  Theoretical  Development 

For  the  system  represented  by  Equation  (2.45),  the  following 

notation  is  used. 

1.  Due  to  the  definition  of  B-j  and  B2  in  Equations  (2.43)  and  (2.44), 
matrices  G-j(s)  and  G2(s)  are  of  the  form 


II  £T) 

II  O 

_ 1 

5  i  =21 

1 

1 

1 

ft  = 

1 

1 

1 

,  «2 

"""T 

1 

§12  i 

2  '  =22 

where  the  partitioning  serves  the  purpose  of  separating  the 
undelayed  and  delayed  control  variables. 

G-j  and  §2  are  n  x  m  matrices  and  if  the  number  of  undelayed 
manipulated  variables  is  £,  the  dimensions  of  the  partitions  are, 
G-j i  =  £X£,  G12  =  (n-Ox£,  G21  =  £x(m-&),  G22  =  (n-a)x(m-a). 

2.  Define  a  matrix  K  as 


0 


0 


where 


K 


mxm  matrix 


i  tin  m  .x  i’.  bn  ,  0 

. 
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K-j  =  (m-£)x£  matrix 
l<2  =  (m-&)x(m-&)  matrix 


3.  The  control  matrix  Gc(s)  is  written  as 


where  the  dimensions  of  the  partitions  are 

Gc  =  mxn ,  Gcl  =  ixz,  (?c2  =  ax(n-Ji),  Gc3  =  (m-a)xJi,  §c4  =  (m-t)x(n-&) 

The  assumption  is  made  that  all  states  are  measurable. 

A  satisfactory  predictor  technique  was  devised,  based  on  the  block 

diagram  in  Figure  2.5  and  by  placing  the  following  two  conditions  on 

matri ces  K  and  Go : 

=  =c 

Condition  1 


Matrix  K  is  selected  such  that 


(2.49) 


(2.50) 


where  I  is  the  (m-Ji)x(m-Ji)  unit  matrix.  This  conditions  is  easily 
satisfied  by  selecting  K-j  =  0,  ]C,  =  J. 

Condition  2 

The  restriction  is  imposed  on  the  partitions  Gc3  and  G^  of 
the  control  matrix  that  they  must  satisfy  the  following  relation: 


=c3  =11  +  =c4  =12 


0 


(2.51) 


It  can  easily  be  shown  by  direct  substitution  that  if  conditions  1  and  2 


■ 


O'. 

5  rid  -oue  >9}39f5  e'r  >t  x 


b-  j 
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are  satisfied,  the  following  equations  hold 


no 

o 

il  o 

n 

no 

(2.52) 

KG  =  G0  G„ 

(2.53) 

=  =c  =2  =c 

K  Gc  =  0 

(2.54) 

By  virtue  of  Equations  (2.52)  -  (2.54)  a  satisfactory  predictor  control 

law  can  be  derived  as  follows. 

Assume  that  a  control  law  is  desired  of  the  form 

u (s )  =  -  Gc  x(s)  -  w(s)  (2.55) 


where 

x(s)  =  measurement  vector  of  dimension  n 
w(s)  =  predictor  output  vector  of  dimension  m 
It  will  be  shown  that  if  w(s)  is  selected  as 

w(s )  =  K  Gc  g(s)  (2-56) 

where 

p(s)  =  predictor  output  of  dimension  n 
then  the  effect  of  the  delayed  control  action  can  be  compensated  by  a 
prediction  with  no  effect  on  the  undelayed  part  of  the  control  action. 
Substitution  of  Equation  (2.56)  into  (2.55)  gives  the  following 
control  law 

u(s)  =  -  Gcx(s)  -  KGcp(s)  (2.57) 

Substitution  of  the  expression  for  u(s)  from  Equation  (2.57)  into 
Equation  (2.47)  gives 


.3Worfol  26  b9vneb  9d  vsd  wsl 
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x(s)  =  -  G1Gcx(s)  -  G^Gggts)  -  G2Gce'asx(s)  -  G2KGce'asp(s ) 

+  GLd(s)  (2.58) 

Combining  Equations  (2.53),  (2.54)  and  (2.58)  gives 

x(s)  =  -  G1Gcx(s)  -  G2Gce'asx(s)  -  G2Gce‘asp(s)  +  GLd(s)  (2.59) 

and  combining  Equations  (2.59),  (2.52),  (2.53)  and  (2.45)  gives 

x (s )  =  -  G1Gcx(s )  +  GLd(s)  -  G2GcGLe'asd(s) 

-  G2G.G2e'2asu(s)  -  G2GLe"asp(s)  (2.60) 

or 

x (s )  =  -  G^s)  +  [In-G2Gce'as]GLd(s) 

-  G2Gce"as[G2e‘asu(s)  +  p(s)]  (2.61) 

Application  of  the  Smith  Predictor  in  the  feedback  loop  around  the 
controller  will  result  in  the  following  expressio n  for  p(s) 

p(s)  =  -  G2  e“asu(s)  +  G2  u(s)  (2£2) 

Combining  Equations  (2.62),  (2.61)  and  utilizing  Equations  (2.52) 

(2.54)  gives 

x(s )  =  -  (G]  +  G2)Gcx(s)  +  [In+G2GcO-e’aS)]G|_d(s)  (2.63) 

or  rearranging 

x(s)  =  [In  +  (G^Ggjy1  [In  +  G2Gc(1  -  e'aS)]GLd(s)  (2.64) 

which  gives  the  closed  loop  expression  for  the  state  vector.  From 


. 


10 
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Equation  (2.64)  it  is  obvious  that  the  characteristic  equation  of  the 
closed  loop  system  is  given  by 

| In  +  %  +  G2)Gc|  =  0  (2.65) 

which  is  also  the  characteristic  equation  of  the  system  without  time 
delays.  It  should  be  noted  that  Equation  (2.65)  can  be  derived  from 
Equation  (2.21)  by  assuming 

Gp  =  £i  +  =2 
-1  +  =2  ”  = 

The  case  of  incomplete  state  measurement  that  is  when 

y (s )  =  C  x(s)  {2.66) 


where 

C  =  rxn  measurement  matrix 

can  also  be  included  in  this  scheme  by  modifying  Equation  (2.52)  as 
follows 

K  G  C  G-j  =  0  (2-67^ 

No  other  modification  of  Equations  (2.52)  -  (2.54)  is  needed  and  the 
resulting  closed  loop  expression  for  x(s)  is 

x(s)  -  [In  +  (G^JScCr1  [In  +  G2GcC(1  -  e-aS)]GLd(s)  (2.68! 

2.4.3  Pi s crete -Time  Systems 

The  discrete  time  case  follows  along  the  same  lines.  Assume 
that  the  state-space  equation  for  the  system  is  given  by 


. 
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x  ( n+1 )  =  *x(n)  +  0-]'u-|(n)  +  ^’u^n-a)  +  =  -(n) 

y(n)  =  x(n) 

where  the  control  vector  is  partitioned  as 


(2.69) 

(2.70) 


(2.71) 


An  equivalent  expression  for  Equation  (2.69)  is  given  by 
x(n+l)  =  $x(n)  +  e,u(n)  +  e2u(n-a)  +  Ad(n) 


(2.72) 


where  e,  and  e2  are  defined  by 


0] 


(2.73) 


§2  =  [0 


—2  ^ 


(2.74) 


and  the  dimensions  are  the  same  as  for  the  continuous  time  case 
Taking  the  Z-transformation  of  Equation  (2.72)  gives 


z  X ( z )  =  ^e1U(z)  +  M02z"aU(z)  +  f;lAD(z) 


(2.73) 


where  M  is  defined  as  in  Equation  (2.36). 

If  we  denote 

Me,  H  e,“  ,  H92  s  e2"  ,  Mi 

then  the  conditions  described  by  Equations  (2.51) 
following  form  in  discrete  time  notation 


A  (2.74) 

(2.54)  take  the 


. 


24 


Figure  2.5  Smith  Predictor  Scheme  for  the  Case  of  Delays  in  Some  Control  Variables 
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£c3lll  +  =c4=12  = 

Hcti"  =  2 
M=Kc  =  Me 
Si"&&c  =  2 


(2.75) 

(2.76) 

(2.77) 

(2.78) 


where  Kc  is  the  discrete  time  control  matrix,  K  is  the  same  matrix  as 
in  Equations  (2.49),  (2.50)  and  the  partitioning  of  matrices  K,  Kc> 
e  ",  e2"  is  done  in  the  same  way  as  in  the  continuous  time  case. 

The  derivation  of  the  closed  loop  expression  for  X(z)  is 
easily  done  by  following  the  similar  continuous  time  procedure  and 

X(z)  is  given  by 


(2.79) 


with  characteristic  equation  given  by 

lz!n  +  %  +  =2  >=Kcl  =  °  (2'8C 

which  is  also  the  characteristic  equation  for  the  system  without  time 
delay. 

2.5  Modification  nf  the  Predictor  Scharr  for  Compensation-a£ 
Inaccurately  Determined  or  Time  Varying  Delays 

2.5.1  General 

An  inaccurate  representation  of  the  delay  in  the  predictor 

loop  can  occur  for  the  following  reasons. 

(1)  Inaccurate  determination  of  the  delay.  . 

(2)  Use  of  an  approximation  (Fade  etc.)  for  the  time  delay  transfer 


function. 
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(3)  Time  varying  delay. 

(4)  Time  delay  which  is  not  an  integer  multiple  of  the  sampling 
interval  in  the  discrete  time  case. 

The  effect  of  approximating  the  actual  delay  by  various  methods 
(Pade,  etc.)  on  the  response  of  the  closed  loop  predictor  system 
has  been  examined  by  Buckley  [9,10]  for  single  variable,  continuous 
time  systems.  In  the  multivariable  discrete  time  case,  a  time  delay 
which  is  not  an  exact  integer  multiple  of  the  sampling  interval  can 
be  handled  by  the  method  used  by  Jacobson  [16]  in  compensating  for 
a  delay  equal  to  half  the  sampling  interval  which  is  a  consequence 
of  the  digital  interface,  [15]. 

2.5.2  Digital  Algorithm  for  the  Determination  of  Inaccurately  Known 
or  Time  Varying  Delays 

In  this  section  an  algorithm  is  derived  to  improve  the 
predictor  performance  when  time  varying  or  inaccurately  determined 
time  delays  occur.  This  digital  algorithm  is  realizable  if  an 
accurate  discrete  time  model  of  the  process  without  delays  is  avail¬ 
able. 

The  derivation  concerns  the  case  of  a  single  delay  in  all 
the  control  variables.  Measurement  delays  are  not  considered  since 
usually  there  is  less  uncertainty  concerning  these  values. 

Analysis  of  the  digital  algorithm 

In  the  discrete  time  Smith  Predictor,  the  time  delay  in  the 
predictor  loop  is  realized  by  storing  the  value  of  the  control  vector 
in  a  table  and  then  selecting  from  this  table  of  past  vectors  the 
appropriate  delayed  control  vector.  In  the  case  of  inaccurate  deter¬ 
mination  of  the  delay,  a  wrong  selection  of  the  control  vector  will 


. 
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adversely  affect  the  control,  especially  in  cases  of  large  sampling 
i ntervals . 

This  method  of  representing  the  time  delay  in  the  predictor 
loop  has  the  advantage  of  making  possible  compensation  for  uncertain 
or  time  varying  delays  according  to  the  following  reasoning. 

Consider  a  system  with  a  delay  in  all  the  manipulated 
variables  equal  to  a  sampling  interval  and  which  can  be  represented 
by  an  accurate,  discrete-time  mathematical  model  as 

x(n)  =  <p  x(n-l)  +  e  u(n-a-l)  +  A  d(n-l)  (2.81) 

where  the  notation  is  the  same  as  for  Equation  (2.22).  For  the  system 
of  Equation  (2.81)  all  state  variables  are  assumed  to  be  measured.  The 
basic  function  of  the  Smith  Predictor  is  the  cancellation  of  the  effect 
of  the  delayed  control  action  and  the  replacement  of  it  by  a  prediction 
of  the  system  response  for  the  system  without  a  time  delay.  In  cases 
of  inaccurate  representation  of  the  delay  in  the  predictor  loop,  the 
cancellation  will  be  incomplete.  However,  if  an  accurate  model  of  the 
undelayed  system  is  available,  the  following  calculations  are  feasible: 

1.  Calculation  of  changes  in  x  due  to  an  unmeasured  load  variable. 

2.  Calculation  of  the  state  vector  for  the  case  where  the  effect  of 
the  delayed  control  action  is  completely  cancelled. 

The  calculation  of  the  state  vector  if  complete  cancellation  of  the 
control  action  is  assumed  can  be  accomplished  by  using  the  model  and 
the  measurement  and  forming  the  following  vector: 


S(n)  =  £  x(n-l )  +  A  d(n-l ) 


(2.82) 
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x(n-l)  =  previous  state  measurement 
d(n-l)  =  load  variable  vector 

The  load  change  can  be  calculated  from  the  first  deviation  of  x  from 
the  steady  state  value.  Any  change  in  the  measurement  can  be  inter¬ 
preted  as  a  load  change,  assuming  that  the  system  is  already  at 
steady  state.  The  measurement  in  perturbed  form  will  then  be  given 
by 


x(n)  =  A  d(n) 

It  should  be  noted  that  only  the  product  A  d(n)  needs  to  be  calculated 
and  that  the  actual  value  of  d(n)  is  not  required.  Under  steady 
state  conditions,  inaccuracies  in  the  predictor  do  not  affect  the 
response  since  the  output  of  the  predictor  will  be  equal  to  zero. 
Combining  Equations  (2.81)  and  (2.82),  the  following  equation  can  be 
wri  tten 

x(n)  -  e  u(n-a-l)  =  S(n)  (2.83) 

According  to  the  scheme  given  in  Figure  2.6,  the  following  procedure  is 
followed: 

(1)  An  approximate  value  for  the  actual  delay  a  is  assumed  and  denoted 
as  a' . 

(2)  The  vector  S(n)  is  calculated  from  Equation  (2.82)  using  the  pre¬ 
vious  measurement,  x(n-l)  and  estimated  or  measured  values  of  the 
disturbance. 

(3)  The  delayed  control  vector  u(n-a'-l)  is  obtained  from  the  stored 
past  values  of  the  control  vector. 
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(4)  The  product  e  u(n-a'-l)  is  calculated  and  subtracted  from  the 
present  measurement  x(n)  and  the  result  compared  to  S(n).  As 
seen  from  Equation  (2.83),  the  result  of  the  forementioned 
subtraction  should  be  equal  to  S(n)  when  a1  =  a. 

If  this  comparison  is  satisfactory,  then  a'  =  a  and  the  vector 
u(n-a'-l)  is  sent  to  the  predictor  loop.  In  the  opposite  case,  the 
search  of  the  table  of  past  values  continues  sequentially  by  assum¬ 
ing  larger  values  of  a'.  In  the  case  of  initially  unknown  but  con¬ 
stant  time  delays,  the  search  can  be  stopped  once  the  value  of  the 
delay  has  been  established.  Assuming  that  the  system  was  initially 
at  steady  state,  determination  of  the  actual  delay  will  be  possible 
only  after  a  time  interval  equal  to  the  actual  time  delay,  since 
until  that  time  the  response  of  the  system  is  open  loop  due  to  the 
time  delay. 

In  the  case  of  slowly  varying  time  delays,  the  previously 
described  search  will  be  continued  but  the  algorithm  has  to  be 
modified  to  include  the  ability  of  distinguishing  between  load 
changes  and  changes  in  the  magnitude  of  the  time  delay.  This  is 
necessary  since  although  the  predictor  control  scheme  does  not 
require  knowledge  of  load  changes,  the  determination  of  the  time 
delay  does  require  knowledge  of  any  load  changes. 

The  algorithm  proceeds  as  follows.  It  is  initially  assumed 
that  any  unexpected  change  in  the  state  vector,  x(n),  is  due  to  a 
change  in  the  value  of  the  time  delay,  a.  Furthermore,  it  is  also 
assumed  that  the  time  delay  is  restricted  to  a  known  range  of  values. 
This  assumption  restricts  the  search  of  the  table  of  past  control 
vectors  to  a  certain  predetermined  range.  If  no  control  vector  is 


■ 
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found  inside  this  range  that  satisfies  the  requirements  of  Equation 
(2.83),  then  a  change  in  the  load  is  assumed  to  have  occurred.  Then, 
using  the  best  available  estimate  of  the  time  delay,  the  respective 
control  vector  is  selected  from  the  table  and  the  new  load  change  is 
calculated  from  the  model.  This  calculation  is  done  by  subtracting 
from  the  current  measurement,  a  calculated  value  of  the  current 
measurement  which  is  formed  as  follows: 

-CAI_(n)  =  +  =  (2.84) 

where  S(n)  is  given  by  Equation  (2.82). 

Subtraction  of  XCA[_(n)  from  the  present  measurement  will  give 

X(n)  -  XCAL(n)  =  A  [d (n -1 )  -  dEST(n-l)]  (2.85) 

which  can  be  solved  for  A  dESJ(n-l)  where  dE$T(n-l)  is  the  estimated 
value  of  the  load  vector.  Since  calculation  of  S(n)  requires  the  value 
of  the  product  A  dEsy(n-l),  n0  calculation  of  the  load  vector  dEgy(n  1) 
is  required  for  the  subsequent  calculations.  This  procedure  is  valid 
as  long  as  a  change  in  the  load  does  not  happen  at  the  same  time  as  a 

change  in  the  time  delay. 

This  algorithm  can  also  be  applied  to  the  case  of  delays  in 
only  some  of  the  manipulated  variables  with  no  serious  modifications 

requi red. 

In  Chapter  3,  simulated  results  are  presented,  that  demon¬ 
strate  the  ability  of  the  algorithm  to  estimate  a  delay  in  the  process 
as  well  as  distinguishing  between  load  changes  and  changes  in  the 
delay,  if  an  accurate  mathematical  model  of  the  system  (without  time 
delays)  is  available.  Simulated  responses  for  a  discrete-time  system 
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Figure  2.6  Predictor  Algorithm  for  Inaccurate  or  Time  Varying  Delay  in  the 
Control  Variables 
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are  also  given  in  Chapter  3,  to  illustrate  the  effect  of  an  inaccurate 
representation  of  the  time  delay  in  the  predictor  loop. 


2.6  A  Numerical  Example  on  the  Multivariable  Predictor  Control 

To  illustrate  the  operation  of  the  multivariable  Smith 
Predictor,  time  delays  were  introduced  into  the  third  order  system 
considered  by  Takahashi  et  al .  [17].  Two  cases  will  be  considered: 

(a)  Time  delay  in  the  control  variables. 

(b)  Time  delays  in  both  the  control  and  output  variables. 

The  system  equations  for  this  example  are: 


x ( t )  =  A  x ( t )  +  B  u(t-a)  +  D  d(t) 


(2.86) 


y(t) 


where : 


C  x(t) 

< 

C-j  x(t)  +  C2  x(t-b) 


for  case  (a) 
for  case  (b) 


(2.87) 
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The  corresponding  discrete  time  model  for  T  =  0.5  can  be 
obtained  in  the  standard  manner,  using  the  analytical  solution  to 
Equation  (2.86)  (Lapidus  and  Luus ,  [18])  and  is  given  in  Equations 

(2.88)  and  (2.89). 
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x(n+l)  =  $  x(n)  +  0  u(n-a)  +  A  d(n) 


(2.88) 


y(n) 


C  x(n) 

< 

C-j  x(n)  +  C2  x(n-b) 


for  case  (a) 
for  case  (b) 


(2.89) 


where 


“0.2840 

0.1310 

0.0606“ 

“0.2760 

0.0148“ 

cj)  - 

0.2620 

0.3440 

0.2620 

0  = 

0.1030 

0.1030 

_0. 0606 
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_0 .01 48 

0. 2760_ 
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A  = 
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_ 0 . 0 1 48 
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A  feedback  control  law  of  the  form 

u(n)  =  -  Kc  y(n) 


(2.90) 


is  desired.  Since  the  Smith  Predictor  algorithm  will  be  used,  Kc  can 
be  designed  using  conventional  design  techniques  for  systems  without 
time  delays,  and  the  state-space  model  in  Equations  (2.86)  and  (2.87) 
with  a  =  b  =  0.  For  example,  a  modification  of  the  direct  synthesis 
method  of  Porter  and  Crossley  [19],  can  be  used  to  design  Kc  so  that  a 
specified  closed-loop  system  matrix,  T,  results. 

This  approach  was  used  to  design  Kc  as 


0.0257  0.0556 

-0.4770  -1.0300 


(2.91) 


corresponding  to  the  closed  loop  system  matrix,  T  =  £  +  ej<cC  where  T 


' 
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was  selected  to  be 


T 


0.2840  0.1310  0.0606 
0.2620  0.2970  0.2000 
0.0606  0  0 


(2.92) 


Transient  responses  of  the  closed-loop  system  including  the 
Smith  Predictor  are  shown  in  Figures  2.7  and  2.8  for  cases  (a)  and  (b), 
respectively.  All  responses  are  for  the  feedback  control  matrix  in 
Equation  (2.91),  the  Smith  Predictor  algorithm  of  Equations  (2.31)  - 
(2.34)  and  a  unit  step  change  in  the  disturbance  vector,  d.  Figure 
2.7  illustrates  the  effect  of  a  delay  in  the  control  variables  on  the 
system  response.  The  open  and  closed-loop  responses  for  the  system 
without  time  delays  are  given  by  curves  5  and  1,  respectively.  When 
the  control  variables  are  delayed  as  in  curves  2-4,  the  system  response 
initially  follows  the  open-loop  response,  then  changes  direction  and 
eventually  approaches  the  closed-loop  response  for  the  undelayed  system. 

Thus,  by  using  the  multivariable  Smith  Predictor  and  a  feed¬ 
back  control  matrix  designed  for  the  system  without  time  delays,  a 
satisfactory  load  response  is  obtained  with  no  offset  between  the 
predictor  response  and  the  response  of  the  undelayed  system. 

In  Figure  2.8  time  delays  in  both  the  control  variables  and 
output  variables  are  considered  for  the  same  disturbance  and  control 
scheme  that  were  used  in  Figure  2.7.  The  responses  are  qualitatively 
similar  to  those  in  Figure  2.7  and  are  judged  to  give  a  satisfactory 
degree  of  control . 
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Figure  2.7  Smith  Predictor  for  System  with  Delay  in  the  Control  Variables 
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Figure  2.8  Smith  Predictor  for  System  with  Delays  in  the  Control  and 
Measurement  Variables 


CHAPTER  THREE 

SIMULATION  RESULTS  FOR  THE  PREDICTOR  SCHEME  APPLIED 
TO  A  DOUBLE  EFFECT  EVAPORATOR 
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3.1  Introduction 

The  multivariable  predictor  algorithm  developed  in  Chapter 
Two  was  simulated  on  an  IBM  1800  Digital  Computer  using  a  5th  order 
state  space  model  of  the  pilot  scale,  double  effect  evaporator  in  the 
Department  of  Chemical  and  Petroleum  Engineering  at  the  University 
of  Alberta.  The  control  algorithm  used  in  the  simulation  was  later 
easily  incorporated  in  the  multivariable  control  program  developed 
by  Newell  [21]  for  experimental  verification  of  the  simulated  results 
on  the  actual  evaporator. 

The  use  of  a  digital  computer  in  the  simulation  and  experi¬ 
mental  studies  of  the  predictor  algorithm  was  ideal  since  the  repre¬ 
sentation  of  the  time  delays  could  be  made  exact  without  approximation. 

3.2  Mathematical  Model 

The  pilot  plant  scale,  double  effect  evaporator  is  represented 

in  the  simplified  schematic  diagram  of  Figure  3.1. 

The  first  effect  is  a  natural  circulation  calandria  type  unit, 
heated  with  a  nominal  2  lb. /min.  of  fresh  steam  and  fed  with  a  nominal 
5  lb. /min.  of  three  percent  triethylene  glycol  by  weight.  First  effect 
vapour  is  used  to  heat  the  second  effect,  an  externally  forced- 
circulation  long  tube  vertical  unit,  which  concentrates  first  effect 
product  to  about  ten  percent.  The  second  effect  is  kept  under  tight 
pressure  control  by  a  vacuum  system  and  condenser.  The  evaporator  has 
been  extensively  modelled  by  Andre  [22],  Newell  [21]  and  Wilson  [20], 
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FIGURE  3.1  Schematic  Diagram  of  Double  Effect  Evaporator  and  Multi  loop  Control  Scheme 
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The  mathematical  model  used  in  the  simulation  and  experimental  runs  is 
due  to  Wilson  [20].  The  nonlinear  model  was  linearized,  with  the 
variables  in  normalized  perturbation  form  and  is  expressed  in  state 
space  notation  by: 

x ( t )  =  A  x(t)  +  B  u(t)  +  D  d(t )  (3.1) 

The  model  used  here  is  a  five  state,  linear  time  invariant 
model  with  matrices  A,  B  and  D  having  the  numerical  values  given  in 
Table  3.1  where  the  state,  disturbance  and  control  vectors  are  denoted 
as 
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where  ,  S',  etc.  denote  normalized  perturbation  variables. 

W,  ’  WUc 

e.g.  W-j  =  - —  (ss  =  steady  state) 
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TABLE  3.1 

Five  State  ,  Continuous  Time  Linearized  Model 
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For  the  digital  algorithm  used  in  the  simulation  and  experimental 
application,  the  state  difference  equation  was  obtained  from  the  con¬ 
tinuous  time  equation  of  Table  3.1  by  using  the  analytical  solution, 

[20],  to  Equation  (3.1)  and  a  time  base  of  64  secs.  The  resulting 
discrete  time  model  is  given  by: 

x(n+l)  =  <j>  x(n)  +  0  u(n)  +  A  d(n)  (3.2) 

with  matrices  <j>,  e  and  A  having  the  numerical  values  given  in  Table  3.2. 
For  the  purpose  of  simulating  measurement  delays,  the  following  output 
equation  was  used 

y (n )  =  C-j  x(n)  +  x(n-b)  (3.3) 

where 

y(n)  =  output  vector  of  dimension  r 
b  =  constant  measurement  delay 
C^,C2  =  constant  matrices  of  dimension  rxn 

The  assumption  was  made  that  all  state  variables  are  available.  For  the 
case  of  delay  in  the  product  concentration,  (C2),  measurement  matrices 
would  be 
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while  for  the  case  of  delay  in  the  control  variables  only, 
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TABLE  3.2 

Discrete  Time  Evaporator  Model  (T  =  64  sec. 
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£t  -  I  •  c2  =  2 

where  I  is  the  5x5  identity  matrix. 

3. 3  Predictor  Algorithm 

For  the  simulated  application  of  the  predictor,  the  equations 
developed  in  Chapter  Two  for  the  Multivariable  Discrete-Time  case  were 

used,  namely.  Equations  (2.22)  -  (2,28). 

The  simulation  of  the  time  delay  in  the  control  variables  was 
accomplished  by  delaying  the  control  action,  and  the  measurement  delay 
by  using  past  values  of  the  stored  state  variable  vector.  The  time 
delay  used  in  the  predictor  algorithm  (Equation  (2.28))  was  achieved  in 
a  similar  fashion  by  forming  a  table  of  past  control  variable  vectors 
and  selecting  the  appropriate  one  according  to  the  value  of  the  time 
delay.  In  all  of  the  simulation  runs  it  was  assumed  that  the  time  delay 
is  an  integer  multiple  of  the  sampling  time.  Throughout  the  simulation 
and  experimental  runs  a  control  law  of  the  form: 

u(n)  =  K(x(n)  +  p(n))  (3-4) 

was  used  where  K  is  a  constant  proportional  control  matrix  and  p(n)  is 
the  predictor  output.  This  control  law  implies  a  change  in  the  sign 
of  the  right  hand  side  of  Equation  (2.25)  and  this  is  the  only  change 
made  to  the  Equations  (2.22)  -  (2.28). 

3.4  Multi  loop  Control 

Multiloop  control  has  been  extensively  applied  to  the  double 
effect  evaporator  at  the  University  of  Alberta,  and  the  theoretical 
considerations  involved  along  with  further  details  can  be  found  in  the 
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theses  by  Newell  [21]  and  Jacobson  [16]. 

A  multi  loop  control  scheme  consists  of  several  single  variable 

control  systems,  i.e.  each  control  variable  is  affected  by  the  changes 
in  only  one  of  the  state  variables.  In  the  evaporator  application  the 
usual  pairing  of  manipulated  and  controlled  variables  is  [16,21,23], 


Control  Variable 

S 
B 


State  Variable 

C2 
W, 


1  1 
B2  W2 

The  control  matrix  for  multiloop  control  was  derived  by  Oliver  [23]  and 


is  given  by 


!5fb 


'  0  0 

3.52  0 

0  0 


0  0 

0  0 

0  15.8 


-4.89 

0 

0 


(3.5) 


A  comparison  of  simulated  open  and  closed  loop  responses  for  the  model 
in  Table  3.1  and  the  multi  loop  control  scheme  of  Equation  (3.5)  is 
presented  in  Figure  3.2,  where  horizontal  arrows  denote  the  initial 
steady  state  values.  As  can  be  seen  in  Figure  3.2,  the  uncontrolled 
(open  loop)  system  is  not  self-regulatory  due  to  the  integrating 
nature  of  the  two  holdups.  The  final  steady  state  values  (i.e.  offsets) 
after  the  +20%  step  change  in  feed  flow  are  given  in  Table  3.3. 

3.4.1  Multi  loop  Control  of  System  With  Time  Delays 

The  actual  double  effect  evaporator  does  not  have  any 
important  time  delays.  Hence,  time  delays  were  arbitrarily  introduced 
into  the  mathematical  model,  and  transient  responses  were  calculated 
for  multi  loop  control  and  different  values  of  the  delays.  In  the 
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TABLE  3.3 

Offsets  Resulting  from  Simulated  Multi  loop  Control  After  a  Step 

Change  of  +20%  in  Feed  Flow 


Vari  able 


W' 


c; 


w: 


Offset  (%) 
5.84 

-0.58 

4.97 

1.52 

-3.40 
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FIGURE  3.2  Simulated  Response  for  Multi  loop  Control  and  a  Step  Change 
of  +20 %  in  Feed  Flowrate 
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simulation  runs  all  states  are  assumed  measurable  and  the  load  change 
consists  of  a  20%  step  increase  in  the  total  feed  flow  starting  at  the 
first  sampling  interval  (i.e.  at  t  =  64  sec.)* 

Results  were  obtained  for  two  cases: 

1)  Time  delay  in  all  the  manipulated  variables,  (Process  Delay). 

2)  A  measurement  delay  in  some  of  the  state  variables. 

In  the  first  case,  a  single  time  delay  was  introduced  in  all  the  mani¬ 
pulated  variables.  Such  a  case  could  arise  in  practice  from  the 
existence  of  long  pneumatic  lines  for  the  actuation  of  pneumatic 
valves  [24].  The  case  of  a  delay  in  some  manipulated  variables  was 
not  examined  in  the  simulated  or  experimental  runs.  However,  a  pre¬ 
dictor  scheme  for  this  case  was  theoretically  derived  and  the  necessary 

conditions  were  presented  in  Chapter  Two. 

The  measurement  delay  case  is  restricted  to  a  time  delay  in 
product  concentration  measurement.  The  behavior  of  the  system  for 
delays  in  the  measurement  of  other  state  variables  is  examined  in  con¬ 
nection  with  multivariable  control  in  §  3.5.  Due  to  the  fact  that  in 
the  double  effect  evaporator  an  on-line  refractometer  is  used  for  con¬ 
tinuous  product  concentration  measurement,  no  significant  delay  exists. 
Therefore,  a  delay  is  introduced  in  the  model  for  purposes  of  study, 
since  delays  in  the  measurement  of  concentration  are  of  common 
occurrence  in  other  processes  and  tend  to  deteriorate  the  control  in 
a  dramatic  fashion.  Examples  of  the  detrimental  effects  of  the  delay 
in  the  concentration  measurement  can  be  found  in  the  study  of  distilla¬ 
tion  column  control  at  the  University  of  Alberta  by  McGinnis  [25]. 

Delay  in  the  Manipulated  Variables 

The  transient  response  of  the  system  was  obtained  under 
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multi  loop  control  for  process  delays  of  64  and  128  sec.  As  shown  in 
Figure  3.3,  the  responses  are  underdamped  for  a  delay  of  64  sec.  and 
remain  underdamped  for  a  delay  of  128  sec.  for  all  states  except  the 
second  effect  holdup  (W2)  which  becomes  unstable.  The  unstable 
behavior  of  W0  is  due  to  the  high  controller  gain  in  the  W2-B2  control 
loop  (i.e.  K34  =  15.8).  Decreasing  this  gain  will  result  in  less 
oscillations  but  the  offset  will  increase.  Furthermore,  the  stable 
closed  loop  response  would  still  tend  to  become  unstable  for  larger 
values  of  the  time  delay.  The  effect  of  lowering  the  controller  gain 
is  considered  in  conjunction  with  measurement  delays  in  Figure  3.8. 
Measurement  Delay 

The  transient  response  of  the  system  for  delays  0,  192  and 
256  sec.  in  the  measurement  of  the  product  concentration  with  multi  - 
loop  control  is  given  in  Figure  3.4.  The  response  remains  under¬ 
damped  for  a  delay  of  192  sec.  with  larger  oscillations  in  C2  than 
in  the  other  state  variables.  For  a  delay  of  256  sec.  the  responses 
become  unstable  for  all  the  state  variables.  Improvement  of  the 
responses  can  again  be  accomplished  by  lowering  the  gains  of  the 
control  matrix  to  a  certain  extent. 

Discussion  of  Results 

A  series  of  simulation  runs  were  conducted  in  order  to  evalu¬ 
ate  the  ability  of  multi  loop  proportional  feedback  control  in  com¬ 
pensating  for  the  detrimental  effects  of  time  delays  in  the  double 

effect  evaporator  model . 

For  both  cases  of  delays  in  the  control  or  measurement  vari¬ 
ables,  the  multi  loop  control  scheme  in  Equation  (3.5)  proved  unsatis¬ 
factory.  For  example,  the  system  is  unstable  for  a  delay  in  all  the 
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FIGURE  3.4  Simulated  Response  for  Multi  loop  Control  and  a  Measurement 
Delay  (+20%  Step  Change  in  Feed  Flow) 
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control  variables  of  128  sec.  or  a  measurement  delay  in  C2  of  192  sec. 
Further  improvement  of  the  response  can  be  accomplished  by  lowering 
the  gains  of  the  feedback  control  matrix  but  this  has  the  well  known 
disadvantage  of  increased  offsets  and  more  sluggish  response. 

3.4.2  Multi  loop  Control  with  a  Predictor 
Introduction 

The  predictor  algorithm,  derived  in  Chapter  Two,  was  simu¬ 
lated  as  discussed  in  §  3.2.2  and  used  together  with  the  multiloop 

control  matrix  given  by  Equation  (3.5). 

Theoretical  considerations  imply  that  the  predictor  response 
will  have  the  same  stability  characteristics  and  follow  the  response 
of  the  undelayed  system  with  the  same  control  matrix.  The  same  two 
cases  are  considered  as  in  the  previous  section. 

Multi  loop  Predictor  with  Delay  in  the  Manipulated  Variables 

The  simulated  response  of  the  evaporator  model  was  obtained 
under  multiloop  proportional  predictor  control  for  delays  of  128,  256 
and  384  sec.  in  all  of  the  manipulated  variables.  The  responses  are 
given  in  Figure  3.5  together  with  the  response  of  the  system  without 
time  delays.  As  it  is  shown  in  Figure  3.5,  the  predictor  response  is 
stable  and  drives  the  response  to  the  one  of  the  undelayed  system. 
Comparison  of  the  responses  in  Figures  3.3  and  3.5  shows  the  improve¬ 
ment  on  the  response  when  the  predictor  is  used,  since  the  system 
without  the  predictor  becomes  unstable  for  a  delay  of  128  sec.  The 
predictor  response  results  in  an  increase  of  the  offset  in  the  first 
effect  holdup  (W-j ) ,  due  to  the  small  gain  of  the  respective  control 
loop  (W-j -B-, ) ,  which  eventually  disappears. 
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Multi  loop  predictor  with  measurement  delay 

The  simulated  response  of  the  evaporator  model  was  obtained 
under  multi  loop  proportional  predictor  control  for  delays  of  128,  384  and 
640  sec.  in  the  measurement  of  the  product  concentration.  The  responses 
are  given  in  Figure  3.6  together  with  the  response  of  the  system 
without  time  delays  and  the  same  control  matrix.  The  stabilizing 
effect  of  the  predictor  is  obvious  by  comparing  the  responses  of 
Figure  3.6  and  Figure  3.4.  The  predictor  response  closely  follows  the 
response  of  the  undelayed  system  with  no  increase  in  the  offset.  The 
responses  of  W-j  and  W2  are  hardly  affected  since  the  delay  is  only  in 
C2  and  multiloop  control  is  used  (i.e.  value  of  S  only  has  a  small 
effect  on  W-j  and  W2).  The  improvement  in  the  response  due  to  the  use 
of  the  predictor  is  more  dramatically  depicted  in  Figure  3.7. 

As  mentioned  in  Section  3.3.1,  the  effect  of  the  time  delays 
can  be  compensated  to  a  certain  extent  by  lowering  the  gains  in  the 
proportional  feedback  control  matrix,  when  a  predictor  is  not  used. 

In  Figure  3.8  a  comparison  is  made  between  the  predictor  response  with 
the  initial  multiloop  proportional  control  matrix  and  the  responses  of 
the  system  without  predictor  but  using  smaller  gain  in  the  C2  -  S 
control  loop  (i.e.  K]5).  The  simulated  responses  correspond  to  the 
evaporator  model  with  a  measurement  delay  of  192  sec.  in  C2-  Since 
the  responses  of  W-, ,  W2  and  C]  are  satisfactory  for  this  case,  the 
gains  of  the  other  control  loops  (W^-B-j  and  W2-B2)  were  left  unchanged. 
As  expected,  the  result  of  the  lower  gain  is  stabilization  of  the 
response.  At  the  same  time,  the  offset  increases  and  the  response 
is  slower  in  reaching  the  steady  state.  Comparison  of  the  predictor 
response  (with  the  initial  control  matrix)  with  the  low  gain  response 
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FIGURE  3.5  Simulated  Multiloop  Predictor  Response  with  Delays  in  the 
Control  Variables  (+20%  Step  Change  in  Feed  Flow) 
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shows  the  advantages  of  the  predictor  scheme.  If  a  larger  delay  is 
used,  these  advantages  will  be  even  more  dramatic. 

3.4.3  Effect  of  Model  Inaccuracy  on  the  Predictor  Response 

The  use  of  an  inaccurate  mathematical  model  for  the  con¬ 
trolled  process  was  examined  for  the  case  of  the  multi  loop  predictor 
control.  The  effect  of  inaccuracies  in  the  mathematical  model  of  the 
double  effect  evaporator  was  examined  by  simulating  an  inaccurate 
delay  in  the  predictor  feedback  loop  around  the  controller.  Due  to 
the  fact  that  the  sampling  interval  used  was  64  sec.  with  delays 
taken  as  integer  multiples  of  this  sampling  time,  the  inaccuracies 
introduced  in  the  simulation  were  severe.  In  the  case  of  delay  in 
the  manipulated  variables,  the  response  was  obtained  for  a  delay  of 
256  sec.  in  all  control  variables  while  the  delay  in  the  predictor 
loop  was  assumed  to  be  320  sec.  As  shown  in  Figure  3.9,  the  response 
is  still  satisfactory  for  C2,  C-j  and  W-j  while  W2  became  unstable. 

This  is  due  primarily  to  the  high  gain  in  the  control  loop  for  W2 
(K34  =  15.8).  Lowering  of  this  gain  from  15.8  to  5.0  resulted  in  a 
stable  response  for  W2  and  the  response  of  the  other  state  variables 
became  even  better.  The  effect  of  using  inaccurate  estimates  of  the 
time  delay  in  the  predictor  is  shown  in  Figure  3.10.  The  low  gain 
(for  W2)  control  matrix  was  used  with  an  actual  delay  of  256  sec.  and 
assumed  delays  in  the  predictor  loop  of  192  and  320  sec.,  respectively 
As  it  is  shown,  the  response  is  worse  for  a  low  estimate  of  the  delay. 

The  case  of  measurement  delay  was  examined  for  a  measure¬ 
ment  delay  of  384  sec.  in  the  product  concentration.  The  initial  feed 
back  control  matrix  in  Equation  (3.5)  was  used  as  well  as  a  predictor 
with  an  estimated  delay  of  320  and  448  sec.,  respectively.  The 
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FIGURE  3.10 


TIME  IN  MINUTES 


Simulated  Multi  loop 
and  a  [Delay  256  sec 
Change  in  Feed  Flow 


Response  With  an  Inaccurate  Predictor 
in  the  Control  Variables  (+20%  Step 
>  K34  =  5.0) 


ii 


60 


CONCN  IN  PC  GLYCOL 


61 


FIGURE  3.12 
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simulated  responses  in  Figure  3.11  indicate  that  the  predictor  response 
remains  essentially  unaffected  for  both  cases  of  over-  and  under¬ 
estimation  of  the  delay  in  the  predictor  loop.  The  effect  of  the 
magnitude  of  the  actual  delay  for  the  same  degree  of  inaccuracy  is 
shown  in  Figure  3.12,  where  the  actual  delays  are  192,  512  and  896 
sec.  while  the  assumed  predictor  delays  are  256,  576  and  960  sec. 

As  can  be  seen  in  Figure  3.12,  the  response  is  far  better  for  small 
values  of  the  actual  delay. 

3.5  Multivariable  Control 

Multivariable  control  of  the  double  effect  evaporator  pilot 
plant  has  been  studied  and  applied  extensively  by  Newell  [21], 

Jacobson  [16],  Wilson  [20],  Hamilton  [26],  and  Oliver  [23].  The  genera¬ 
tion  of  an  appropriate  proportional  feedback  control  matrix  is  done  by 
applying  discrete  dynamic  programming  techniques  to  a  formulation  based 
on  a  quadratic  performance  index,  and  the  standard  linear,  time 
invariant,  state  space  model.  The  model  used  in  the  simulation  of 
multivariable  control  is  the  discrete  fifth  order  state  space  model 
given  in  Section  3.2.  The  proportional  feedback  control  matrix  used 
is  due  to  Wilson  [20]  and  is  derived  for  heavy  weighting  of  the  product 
concentration  in  the  performance  index,  it  is  given  in  Table  3.4. 

All  state  variables  are  assumed  to  be  measured  as  in  the 
multi  loop  case.  Examination  of  this  proportional  feedback  control 
matrix  indicates  that  steam  flowrate  (S)  depends  mainly  on  product 
concentration,  (C2),  and  the  product  flowrate  (B2)  depends  mainly  on 

and  the  second  effect  liquid  holdup,  W ^ . 

The  simulated  closed  loop  response  for  this  control  configura- 

tion  is  given  in  Figure  3.13  together  with  the  open  loop  response  for  a 
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TABLE  3.4 

Multivariable  Control  Scheme  and  Offsets 
Multivariable  proportional  feedback  control  matrix 
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Steady  state  offsets  for  step  change  of  +20%  in  feed  flow 
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FIGURE  3.13  Simulated  Multivariable  Response  for  a  +20%  Step  Change 
in  Feed  Flowrate 
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step  change  in  feed  flowrate  of  +20%,  introduced  on  the  second  sampling 
interval.  The  final  steady  state  values  of  the  perturbation  variables 

are  also  given  in  Table  3.4. 

3.5.1  Multivariable  Control  of  System  with  Time  Delays 

As  mentioned  in  Section  3.4.1,  the  double  effect  evaporator 
does  not  have  any  important  time  delays.  For  the  simulation  runs  under 
multivariable  control  process,  measurement  and  combination  of  process 
and  measurement  delays  were  assumed  and  the  load  responses  were  for  a 
step  change  of  +20%  in  the  feed  flow,  introduced  on  the  second  sampling 

interval . 

Delay  in  the  Manipulated  Variables 

The  simulated  response  of  the  double  effect  evaporator  is 

given  in  Figure  3.14  for  delays  of  0,  64  and  128  sec.  in  all  the  mani - 
pulated  variables.  As  can  be  seen  from  Figure  3.14,  the  system  becomes 
unstable  for  a  delay  of  128  sec.  Comparison  of  Figures  3.14  and  3.3 
shows  that  the  multiloop  control  scheme  gives  a  more  stable  response 
due  to  the  lower  controller  gains. 

Measurement  Delay 

The  simulated  response  of  the  system  in  the  presence  of 
measurement  delays  is  given  in  Figure  3.15.  A  delay  of  64,  128  and  192 
sec.  in  the  measurement  of  the  product  concentration,  (0,),  was  con¬ 
sidered.  The  system  response  gradually  becomes  unstable  with  a  delay 
of  128  sec.  and  is  obviously  unstable  for  a  delay  of  192  sec.  As  in 
the  case  of  delays  in  the  manipulated  variables,  the  response  of  the 
system  with  a  delay  is  better  for  multiloop  control  than  it  is  for 
multivariable  control  (cf.  Figures  3.15  and  3.4).  In  Figure  3.16  the 
simulated  response  is  given  for  the  case  of  a  delay  of  192  sec.  in  the 
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FIGURE  3.14  Simulated  Multivariable  Response  with  a  Delay  in  the  Control 
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CONCN  IN  PC  GLYCOL 


67 


0 


30  GO  30 


.1  SO 


TIME  IN  Ml NUT EG 


FIGURE  3.15  Simulated  Multivariable  Response  With  Measurement  Delay  in 
C2  (+20%  Step  Change  in  Feed  Flow) 
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FIGURE  3.16 
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Simulated  Multivariable  Response  with  Measurement  Delay 
192  sec.  in  W-,  and  W2  (+20%  Step  Change  in  Feed  Flow) 
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measurement  of  W-j.  The  overall  response  is  better  than  for  the  case  of 
delay  of  the  same  magnitude  in  concentration  measurement  due  to  the 
larger  weighting  of  C2  in  the  derivation  of  the  optimal  feedback  control 
matrix.  As  is  shown  in  Figure  3.16,  the  measurement  delay  in  the  first 
effect  holdup  (W-j)  affects  the  concentration  more  than  a  delay  in  the 
measurement  of  W2,  as  would  be  expected  from  physical  intuition  (i.e. 

B-|  affects  C2  while  B2  does  not). 

3.5.2  Multivariable  Predictor  Control 

The  same  predictor  algorithm  as  in  the  multiloop  case  is  used, 
the  only  difference  being  the  use  of  the  proportional  feedback  control 
matrix  given  in  Table  3.4.  Three  cases  are  considered:  delay  in  the 
manipulated  variables,  measurement  delays  and  a  combination  of  both 
types  of  delays. 

Multivariable  Predictor  with  Delays  in  the  Manipulated  Variables 

The  simulated  response  of  the  evaporator  model  under  multi - 
variable  predictor  control  was  obtained  for  process  delays  of  64,  128 
and  256  sec.  As  is  shown  in  Figure  3.17,  the  response  is  stable,  and 
eventually  approaches  the  response  of  the  undelayed  system  except  for 
the  first  effect  holdup  (W-j)  that  shows  the  same  apparent  offset  as  was 
observed  in  the  multi  loop  simulated  runs. 

Multivariable  Predictor  with  Delays  in  the  Measurements  and/or 
Manipulated  Variables 

The  simulated  response  for  the  evaporator  with  multivariable 
predictor  control  and  delays  of  128,  384  and  640  sec.  in  the  measure¬ 
ment  of  the  product  concentration,  (C2),  is  given  in  Figure  3.18 
together  with  the  response  of  the  system  without  delays.  For  predictor 
control  the  response  is  stable  and  follows  closely  the  response  of  the 
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FIGURE  3.18  Simulated  Multivariable  Predictor  Response  with  Measurement 
Delay  in  C2  (+20%  Step  Change  in  Feed  Flow) 


CGNCN  IN  PC  GLYCOL 


72 


TIME  IN  MINUTES 

FIGURE  3.19  Simulated  Multivariable  Predictor  Response  with  a  Delay  of 
128  sec  in  the  Control  Variables  and  a  Measurement  Delay 
in  C2  (+20%  Step  Change  in  Feed  Flow) 
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system  without  delays.  No  difference  in  offsets  is  obtained  between 
the  predictor  response  and  the  response  of  the  undelayed  system.  In 
Figure  3.19  the  response  is  given  for  delays  in  both  manipulated  and 
measured  variables.  The  introduction  of  the  delay  in  the  manipulated 
variables  resulted  in  increased  offset  in  W-|  similar  to  that  obtained 
for  delay  in  only  the  manipulated  variables,  (Figure  3.5);  otherwise, 
the  response  is  similar  to  the  one  given  in  Figure  3.18  for  measure¬ 
ment  delays. 

Effect  of  Model  Inaccuracy  on  the  Predictor  Response 

The  effect  of  model  inaccuracy  was  simulated  (as  in  the 
multi  loop  case)  by  using  inaccurate  representation  of  the  time  delay 
in  the  predictor  loop.  The  case  of  delayed  control  variables  proved 
to  be  extremely  sensitive  to  inaccurate  estimates  of  the  delay.  The 
simulated  response  is  given  in  Figure  3.20  for  an  actual  delay  of 
512  sec.  with  inaccuracy  of  t  64  sec.  in  the  assumed  predictor  loop 
delay.  The  response  is  unstable  and  is  generally  worse  for  a  low 
estimate  (448  sec.)  than  a  high  estimate  (576  sec.)  of  the  delay. 

Similarly  for  a  measurement  delay  in  C ^ »  The  effect  of  an 
inaccuracy  of  t  64  sec.  in  the  assumed  predictor  loop  delay  results 
in  unstable  response  (not  shown).  Improvement  of  the  response  is 
obtained  for  larger  vslues  of  the  time  delay  and  the  same  degree  of 

inaccuracy  as  shown  in  Figure  3.21. 

It  should  be  noted  that  conclusions  drawn  from  the  simu¬ 
lated  responses  of  Figures  3.20  and  3.21  should  not  be  generalized, 
since  the  dependence  on  sampling  time  and  delay  is  an  important  factor; 
eg.  faster  sampling  time  will  result  in  less  sensitivity  for  the  same 
inaccuracy  of  t  1  sampling  interval  in  the  predictor  loop  delay, 
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assuming  that  the  delay  is  taken  as  an  integer  multiple  of  the  sampling 
time.  For  instance,  the  sensitivity  will  be  different  for  an  error  of 
+  50%  in  the  representation  of  the  delay,  in  the  predictor  loop,  accord¬ 
ing  to  the  sampling  time  used.  As  a  general  conclusion,  though,  it  can 
be  stated  that  in  the  evaporator  application  the  multivariable  predictor 
is  more  sensitive  to  inaccuracies  in  the  predictor  loop  than  the  multi- 
loop  predictor  due  to  the  large  controller  gains  in  multivariable 
control  (cf.  Equation  (3.5)  and  Table  3.4). 

3.6  The  Effect  of  an  Inaccurate  Time  Delay  in  the  Predictor  Loop  and__a 
Compensating  Algorithm 

In  this  section  the  compensating  algorithm  for  inaccurate  or 
time  varying  delays  (see  Section  2.5)  is  applied  to  the  double  effect 
evaporator  model  under  multiloop  predictor  control.  A  delay  of  192  sec. 
is  introduced  in  all  the  manipulated  variables  while  the  predictor  is 
initially  using  a  delay  of  128  sec.  According  to  the  theory  developed 
in  §  2.5,  the  compensating  algorithm  will  start  a  search  of  the  stored 
past  control  values  according  to  the  algorithm  developed  in  §  2.5.  The 
search  of  the  past  control  values  starts  from  the  initial  assumption  of 
a  delay  of  128  sec.  and  continues  for  larger  values  of  the  delay.  To 
demonstrate  the  effect  of  an  unmeasured  change  in  the  load  variable,  the 
search  of  the  past  control  values  continues  at  every  sampling  interval, 
although  the  delay  is  assumed  constant,  and  the  initial  disturbance  of 
+20%  in  feed  flow  is  changed  after  80  sampling  intervals  to  +10%. 

The  simulated  response  is  given  in  Figure  3.22  together  with 
the  responses  for  the  ideal  case  where  the  predictor  used  the  actual 
value  of  the  time  delay,  (i.e.  192  sec.)  and  is  aware  of  the  feed  flow- 
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responses  for  the  ideal  predictor  case.  The  actual  value  of  the  delay 
is  determined  after  a  time  interval  approximately  equal  to  the  time 
delay.  The  predictor  has  knowledge  of  the  initial  disturbance  of  +20 % 
but  not  of  the  load  change  from  +20%  to  +10%.  In  the  example  of  Figure 
3.22,  since  the  actual  delay  is  constant,  the  search  should  stop  as 
soon  as  the  exact  delay  was  found  and  this  value  used  in  the  predictor 
loop.  The  predictor  then  would  function  according  to  the  original 
predictor  scheme  without  further  use  of  the  compensating  algorithm. 
However,  the  search  was  continued  at  each  sampling  instant  in  order  to 
demonstrate  the  ability  of  the  algorithm  to  sense  the  change  in  the 
disturbance  according  to  the  logic  developed  in  §  2.5. 

It  should  be  noted  that  under  steady  state  conditions  any 
value  of  the  delay  in  the  predictor  loop  will  not  affect  the  response. 
Also,  the  continuous  search  is  necessary  only  if  time  varying  delays 

occur. 

3.7  Discussion  of  Simulation  Results 

The  effect  of  the  time  delays  on  the  simulated  response  of  the 

double  effect  evaporator  was  proved  to  result  in  very  bad  response  for 
both  conventional  multiloop  and  multivariable  feedback  control  schemes. 
The  system  becomes  unstable  for  small  values  of  the  time  delay.  Com¬ 
pensation  for  the  time  delay  effect  can  be  accomplished  by  lowering  the 
controller  gains  but,  inevitably,  results  in  larger  offsets  and  slower 
responses.  This  procedure  does  not  always  result  in  a  satisfactory 
response  especially  when  large  delays  are  involved.  The  response  of  the 
delayed  system  was  generally  better  with  multiloop  control  than  with 
multivariable  control  for  the  control  matrices  used  in  the  simulation 
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Introduction  of  the  Smith  Predictor  into  the  control  scheme 
significantly  improved  the  response  for  delays  in  the  measurement 
and/or  the  control  variables.  Two  control  schemes  were  examined  in 
the  Smith  Predictor  case,  multi  loop  and  multivariable  predictor 
control.  In  both  cases,  the  detrimental  effects  of  the  time  delay 
were  eliminated,  even  for  very  large  delays.  The  response  of  the 
system  follows  the  response  of  the  undelayed  system  with  a  small 
increase  in  the  offset  of  the  first  effect  level  in  the  case  of  delays 
in  the  manipulated  variables.  This  increase  in  the  offset  is  due  to 
the  small  controller  gains  and  eventually  decreases  to  the  respective 
value  for  the  undelayed  system.  The  Smith  Predictor  does  not  intro¬ 
duce  any  increase  in  offset  and  this  can  be  seen  in  the  numerical 

example  in  §  2.6. 

Inaccuracies  in  the  process  model,  expressed  as  an  inaccurate 
representation  of  the  time  delay  in  the  predictor  loop  badly  affected 
the  response  when  multivariable  predictor  control  was  used,  while  the 
multiloop  predictor  scheme  demonstrated  less  sensitivity  to  modelling 
errors.  Finally,  use  of  smaller  controller  gains  improves  the  pre- 
dictor  response  in  the  case  of  large  modelling  errors. 
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CHAPTER  FOUR 

THE  MULTIVARIABLE  SMITH  PREDICTOR  APPLIED  TO  A  DOUBLE 
EFFECT  EVAPORATOR:  EXPERIMENTAL  RESULTS 

4.1  Introduction 

In  this  chapter  an  experimental  application  of  the  multi- 
variable  Smith  Predictor  algorithm  is  presented.  The  experimental 
system  consisted  of  a  pilot  scale  double  effect  evaporator,  inter¬ 
faced  to  an  IBM  1800  data  acquisition  and  control  computer.  The 
theoretical  development  of  the  predictor  algorithm  as  well  as  simulated 
results  have  been  presented  in  the  previous  two  chapters.  This  chapter 
is  concerned  with  the  implementation  of  the  predictor  algorithm  on  the 
experimental  system  and  with  the  results  obtained. 

4.2  Control  System 

The  implementation  of  the  predictor  algorithm  was  achieved 
through  a  modification  of  the  multivariable  control  program  developed 
by  Newell  [21].  The  required  steps  for  the  implementation  of  multi - 
variable  control  are  given  in  the  block  diagram  in  Figure  4.1.  An 
IBM  1800  data  acquisition  and  control  computer  is  interfaced  to  a 
pilot  plant  double  effect  evaporator,  shown  diagramatically  in  Figure 
4.2,  through  analog  to  digital  converters,  digital  to  analog  converters, 
and  required  transducers.  Detailed  descriptions  of  the  equipment  as 
well  as  the  required  electronic  instruments  etc.  can  be  found  in  the 

theses  of  Newell  [21],  Fehr  [27]  and  Andre  [22]. 

The  multivariable  control  program  developed  by  Newell  utilizes 

the  standard  Direct  Digital  Control  (DDC)  package  [28]  for  assessing 
measurements  of  the  state  variables  and  implementing  control  through 
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FIGURE  4.1  Schematic  Diagram 
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FIGURE  4.2  Schematic  Diagram  of  Double  Effect  Evaporator  and  Multiloop  Control  Scheme 
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adjustments  of  setpoints  in  the  DDC  loops. 

The  standard  multivariable  control  program  shown  in  block 
diagram  in  Figure  4.1  performs  model  calculations  for  estimating  the 
unmeasured  first  effect  concentration  (C^ )  and  performs  all  the 
necessary  conversions  of  analog  signals  to  engineering  units  and 
vice  versa.  The  control  program  which  is  written  in  FORTRAN  except 
for  assembler  subroutines  used  for  communications  with  DDC,  is  trans¬ 
ferred  from  disk  storage  to  core  for  execution  at  specified  time 
intervals  designated  by  a  high  priority  timer  interrupt.  In  the 
experimental  runs  presented  in  this  chapter,  the  control  interval 
was  64  seconds.  Although  the  standard  multivariable  program  requires 
from  4  to  7  seconds  execution  time,  the  modifications  introduced  for 
the  predictor  scheme  considerably  increased  the  size  of  the  computer 
program.  This  resulted  in  an  increase  of  the  execution  time  since 
core  storage  limitations  necessitated  an  increase  in  disk  swapping. 

The  modifications  for  the  predictor  application  consisted 

of  the  following  parts: 

1.  The  evaporator  pilot  plant  does  not  have  any  significant  time 
delays  and  no  provision  for  such  a  case  was  included  in  the  multi- 
variable  program.  The  time  delay  was  easily  simulated  in  the 
multivariable  program  by  disk  storage,  and  subsequent  retrieval, 
of  the  measured  states.  This  procedure  refers  to  simulation  of 
measurement  delays;  the  model  calculations  for  the  first  effect 
concentration  proceeded  as  in  the  standard  multivariable  control 

program. 

2.  The  program  included  the  discrete-time  Smith  Predictor  algorithm 
that  was  used  in  the  simulation  study.  This  was  easily  incorporated 
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since  the  standard  multivariable  program  stores  values  of  the  past 
control  vectors.  Consequently,  these  values  could  be  retrieved  from 
disk,  converted  to  perturbation  form,  and  inserted  in  the  calculations 
of  the  predictor  algorithm. 

The  case  of  delays  in  the  manipulated  variables  can  also 
be  easily  included  in  the  multivariable  control  program  by  sending  a 
delayed  control  signal  to  the  process. 

4.3  Experimental  Runs 

The  predictor  control  scheme  was  experimentally  applied  for 
the  case  of  a  measurement  delay  in  product  concentration.  All  experi¬ 
mental  runs  utilized  a  step  change  of  +20%  in  feed  flow  at  the 
initiation  of  the  run.  The  initial  steady  states  were  between  9  and 
10%  in  product  concentration.  Two  values  for  the  measurement  delay 
were  used,  256  and  512  sec.  Due  to  equipment  difficulties,  multi - 

variable  control  runs  were  not  made. 

The  behavior  of  the  undelayed  system  under  multiloop  control 

is  shown  in  Figure  4.3.  The  control  matrix  used  in  all  the  multiloop 
experimental  runs  is  the  one  given  in  Table  3.1  for  the  simulation 
studies.  The  initial  and  final  steady  state  values  of  the  four 
measured  states  are  shown  in  Table  4.1.  As  indicated  in  Figure  4.3, 
the  +20%  step  change  in  the  feed  flow  results  in  significant  offsets 
in  W1  and  C2  since  only  proportional  feedback  control  is  being  used. 
In  Figure  4.4  the  response  is  given  for  multiloop  predictor  control 
with  a  time  delay  of  256  sec.  in  the  measurement  of  product  concentra 
tion.  This  delay  was  introduced  as  described  in  Section  4.2  and  the 
same  control  matrix  was  used  as  for  the  undelayed  case  of  Figure  4.3. 
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Comparison  of  the  state  and  the  control  variables  in  Figures  4.3  and 
4.4  shows  that  the  predictor  response  is  more  oscillatory  but  follows 
the  response  of  the  undelayed  system  quite  well  with  no  increase  in 
the  offset.  This  is  more  clearly  shown  in  Figure  4.5  where  an  overplot 
is  given  of  the  experimental  responses  corresponding  to  Figures  4.3 
and  4.4.  In  Figure  4.6  the  simulated  and  experimental  multiloop 
predictor  responses  are  compared  for  a  measurement  delay  of  256  sec. 
in  C2.  Figure  4.6  indicates  that  the  two  responses  agree  quite  clearly. 
In  Figure  4.7  the  experimental  multiloop  predictor  response  is  given 
for  a  measurement  delay  of  512  sec.  in  The  response  of  C2  is  again 
satisfactory  and  only  W2  shows  an  increase  in  the  degree  of  oscillations. 
In  Figure  4.8  the  simulated  and  experimental  multi  loop  predictor 
responses  are  compared  for  a  measurement  delay  of  516  sec.  in  C2. 

The  experimental  responses  for  B2  shown  in  Figures  4.3,  4.4 
and  4.7  are  somewhat  noisier  than  in  previous  studies  [21,23].  This 
was  due  to  a  noisy  flow  transmitter  which  was  corrected  later. 

4.4  Discussion  of  Results 

The  experimental  results  confirm  the  conclusions  drawn  from 
the  simulation  studies.  The  predictor  performed  satisfactorily  for 
delays  in  the  measurement  of  product  concentration  as  large  as  512  sec. 
The  predictor  response  follows  closely  the  response  of  the  system 
without  delays  without  undue  oscillations  or  increase  in  offsets. 
Furthermore,  the  experimental  runs  demonstrate  one  of  the  main  advan¬ 
tages  of  the  Smith  Predictor  method,  namely,  that  measurement  of  the 
disturbance  is  not  required. 

Another  conclusion  drawn  from  the  experimental  runs  is  that 
modelling  errors  that  inevitably  exist,  did  not  significantly  affect 
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the  predictor  response  and  this,  again,  is  in  agreement  with  the  simu¬ 
lation  studies. 

Finally,  physical  constraints  on  the  control  variables  which 
were  not  incorporated  in  the  simulation  studies  did  not  have  adverse 
effects  on  the  experimental  predictor  response. 
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FIGURE  4.6 


mpari son  of  Experimental  and  Simulated  Multi  loop  Predictor 
sponses  for  Measurement  Delay  of  256  sec.  in  C2  (  2U/0 
ep  Change  in  Feed  Flow) 
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FIGURE  4.7  Experimental  Multiloop  Predictor  Response  wi 
(+20%  Step  Change  in  Feed  Flow) 
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FIGURE  4.8  Comparison  of  Experimental  and  Simulated  Multi loop  Predictor 
Responses  for  Measurement  Delay  of  516  sec.  in  (+20/, 

Step  Change  in  Feed  Flow) 


Table  4.1 

Steady  States  for  Experimental  Multi  loop  Runs 
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CHAPTER  FIVE 
CONCLUSIONS 

The  classical  Smith  Predictor  Method  (Smith,  [5,6])  for  single 
variable  control  systems  has  been  extended  to  a  class  of  linear  multi  - 
variable  systems.  Multivariable  Smith  Predictors  have  been  derived  for 
both  continuous -time  and  discrete-time  systems  which  contain  time 
delays  in  either  the  control  variables  and/or  the  output  variables. 

An  important  advantage  of  the  multivariable  Smith  Predictor  is  that  it 
eliminates  time  delays  from  the  characteristic  equation  of  the  closed 
loop  system.  This  allows  the  control  system  designer  to  choose  from 
the  wide  variety  of  synthesis  techniques  available  for  systems  without 
time  delays  as  opposed  to  the  much  smaller  number  of  techniques  which 
are  applicable  to  systems  containing  time  delays.  Another  advantage  of 
the  multivariable  Smith  Predictor  is  that  the  discrete-time  algorithm 
in  Equations  (2.25)  -  (2.28)  can  be  easily  implemented  on  a  real-time 
digital  computer. 

Simulation  and  experimental  runs  were  carried  out  to  investi¬ 
gate  the  performance  of  the  predictor  algorithm  on  a  double  effect 
evaporator  pilot  plant.  The  simulation  runs  considered  multiloop  and 
multivariable  control  scheme  while  the  experimental  runs  with  the 
evaporator/ IBM  1800  control  configuration  confirmed  the  simulation 
results  for  multi  loop  control.  Simulation  runs  also  examined  the  effect 
of  modelling  errors  on  the  predictor  response  with  multi  loop  and 
multivariable  control. 

Furthermore,  systems  with  time  delays  in  some  of  the  control 
variables  and  time  varying  or  inaccurately  determined  delays  were 
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theoretically  investigated  and  algorithms  were  developed  for  handling 
these  more  complex  systems. 
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NOMENCLATURE 


Alphabetic 

a  time  delay  in  control  variables 

state  coefficient  matrix 
time  delay  in  output  variables 
control  coefficient  matrix 
first  effect  bottoms  flow 
second  effect  bottoms  flow 
output  coefficient  matrix 
C1  first  effect  concentration 

C2  second  effect  concentration 

Cp  feed  concentration 

d  load  variable 

d  load  vector 

D  z-transform  load  vector 

D  disturbance  coefficient  matrix 

e  error 

e  error  vector 


F  feed  flow 

G  transfer  function 

G  transfer  function  matrix 

h1  first  effect  enthalpy 

hF  feed  enthalpy 

H  output  transfer  function 

I  integral  control 

I  identity  matrix 

k.  .  elements  of  control  matrix 

U 


' 


/  r 


II  ^  II 


Nomenclature  (continued) 


1 

m 

M 

n 

P 

P 

P 

PI 

PID 

r 

r 

s 

S 

S 

SI 

t 

T 

u 

u 

u 

w 


control  matrix 

constant  matrix  in  Figure  2.5 

undelayed  control  vector  dimension 

control  vector  dimension 

matrix  defined  in  Equation  (2.36) 

state  vector  dimension 

load  vector  dimension 

predictor  output  vector 

z-transform  of  predictor  output  vector 

proportional -integral  control 

proportional-integral -derivative  control 

output  vector  dimension 

setpoint  vector 

Laplace  operator 

steam  flow 

vector  defined  in  Equation  (2.82) 

sampling  interval 

time 

discrete  state  coefficient  matrix 
control  variable 
control  vector 

z-transform  of  control  vector 
vector  defined  in  Equation  (2.55) 
first  effect  holdup 
second  effect  holdup 
state  variable 


x 
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Nomenclature  (continued) 


x  state  vector 

X  z-transform  of  state  vector 

y  output  variable 

y  output  vector 

Y  z-transform  of  output  vector 

z  z-transform  operator 

Greek 

discrete  disturbance  coefficient  matri 
discrete  control  coefficient  matrix 
discrete  state  coefficient  matrix 


Subscripts 

c  controller 

cj  partitions  of  control  matrix 

CAL  calculated 

EST  estimated 

FB  feedback 

L  load 

p  process 

ss  steady  state 

Superscripts 

1  perturbation  variable 

time  derivative 
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